Pion Interferometry in Au+Au and Cu+Cu Collisions at a/snn = 62.4 and 200 GeV 
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We present a systematic analysis of two-pion interferometry in Au+Au collisions at -^snn = 62.4 
GeV and Cu+Cu collisions at ^^s^ = 62.4 and 200 GeV using the STAR detector at RHIC. The 
multiplicity and transverse momentum dependences of the extracted correlation lengths (radii) are 
studied. The scaling with charged particle multiplicity of the apparent system volume at final inter- 
action is studied for the RHIC energy domain. The multiplicity scaling of the measured correlation 
radii is found to be independent of colliding system and collision energy. 

I. INTRODUCTION or density [H strongly interacting matter will be in a 



One of the definitive predictions of quantum chromo- 
dynamics (QCD) is that at sufficiently high temperature 
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state with colored degrees of freedoni, i.e. quarks and 
gluons. The central goal of the experiments with rela- 
tivistic heavy ion collisions is to create and study this hy- 
pothesized form of matter, called the quark-gluon plasma 
(QGP), which might have existed in the microsecond old 
universe. Numerous experimental observables have been 
proposed as signatures of QGP creation in heavy ion col- 
lisions [2] ■ One of these predictions is based on the expec- 
tation that the increased number of degrees of freedom 
associated with the color deconfined state increases the 
entropy of the system which should survive subsequent 
hadronization and freeze-out (final interactions). The in- 
creased entropy is expected to lead to an increased spatial 
extent and duration of particle emission, thus providing 
a significant probe for the QGP phase transition 0, 0] . 

The information about the space-time structure of the 
emitting source can be extracted with intensity interfer- 
ometry techniques This method, popularly known 
as Hanbury Brown and Twiss (HBT) correlations, was 
originally developed to measure angular sizes of stars [1] . 
The momentum correlations of the produced particles 
from hadronic sources however include dynamical as well 
as interference effects, hence the term femtoscopy is 
more appropriate. The primary goal of femtoscopy, per- 
formed at mid-rapidity and low transverse momentum, 
is to study the space-time size of the emitting source and 
freeze-out processes of the dynamically evolving collision 
system. Femtoscopic correlations have been successfully 
studied in most of the heavy ion experiments (see Q for 
a recent review). 

Experimentally, the two-particle correlation function is 
the ratio, 



C{q,K) = 



A{q,K) 
B{q,K) 



(1) 



where A{q, K) is the distribution of pairs of particles with 
relative momentum q = pi ^ P2 and average momentum 
K = {pi +p2)/2 from the same event, and B{q,K) is 
the corresponding distribution for pairs of particles taken 
from different events 0, ■ The correlation function is 
normalized to unity at large q. With the availability of 
high statistics data and development of new techniques, 
it has become possible to measure three-dimensional de- 
compositions of q [111, m, [l^l , providing better insight 
into the collision geometry. 

Previous femtoscopic measurements at RHIC in 
Au-I-Au coUisions at ^/iiS" = 130 GeV [H, [11] and 
200 GeV [H, [l3 obtained qualitatively similar source 
sizes. However, detailed comparisons with smaller col- 
liding systems and energies are required in order to un- 
derstand the dynamics of the source during freeze-out. 
The crucial information provided from such femtoscopic 
studies with pions will help to improve our understanding 
of the reaction mechanisms and to constrain theoretical 
models of heavy ion collisions [H, [H HO, [HI, [H [H, H, 

In this paper we present a systematic analysis of two- 
pion interferometry in Au-t-Au collisions at ^/snn — 



62.4 GeV and Cu-|-Cu collisions at ^^snn = 62.4 GeV and 
200 GeV using the Solenoidal Tracker at RHIC (STAR) 
detector at the Relativistic Heavy Ion Collider (RHIC). 
The article is organized as follows : Section |TT] explains 
the detector set-up, along with the necessary event, par- 
ticle and pair cuts. In Section [IIII the analysis and con- 
struction of the correlation function is discussed. The 
presented results are compared with previous STAR mea- 
surements for Au-|-Au collisions at y'SNN = 200 GeV in 
Section IIVI This section also includes a compilation of 
freeze-out volume estimates for all available heavy ion 
results from AGS, SPS and RHIC. Section IVl contains a 
summary and conclusions. 



II. EXPERIMENTAL SETUP, EVENT AND 
PARTICLE SELECTION 

A. The STAR detector and Trigger details 

The STAR detector [26*1, which has a large acceptance 
and is azimuthally symmetric, consists of several detec- 
tor sub-systems and a solenoidal magnet. In the present 
study the central Time Projection Chamber (TPC) [l^l 
provided the main information used for track reconstruc- 
tion. It is 4.2 m long and 4 m in diameter. The TPC 
covers the pseudo-rapidity region \r]\ < 1.8 with full az- 
imuthal coverage (-tt <4><tt ). It is a gas chamber filled 
with PIO gas (10% methane, 90% argon) with inner and 
outer radii of 50 and 200 cm, respectively, in a uniform 
electric field of ^ 135 V/cm. The paths of the particles 
passing through the gas are reconstructed from the re- 
lease of secondary electrons that drift to the readout end 
caps at both ends of the chamber. The readout system is 
based on multi-wire proportional chambers with cathode 
pads. There are 45 pad-rows between the inner and outer 
radii of the TPC. 

A minimum bias trigger is obtained using the charged 
particle hits from an array of scintillator slats arranged in 
a barrel, called the Central Trigger Barrel, surrounding 
the TPC, two Zero-Degree Calorimeters (ZDCs) [28j] at 
±18 m from the detector center along the beam line, and 
two Beam-Beam Counters. The ZDCs measure neutrons 
at beam rapidity which originate from the break-up of 
the colliding nuclei. The centrality determination which 
is used in this analysis is the uncorrected multiplicity of 
charged particles in the pseudo-rapidity region |?7| < 0.5 
i^Ih"^) as measured by the TPC. 



B. Event and Centrality Selection 

For this analysis we selected events with a collision 
vertex within ±30 cm measured along the beam axis from 
the center of the TPC. This event selection is applied to 
all the data sets discussed here. 

The events are further binned according to collision 
centrality which is determined by the measured charged 
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TABLE I: Collision centrality selection in terms of percentage 
of total Au-Au inelastic cross-section, number tracks in TPC, 
average number of participating nucleons and average number 
of binary nucleon-nucleon collisions for Au+Au at ^snn ~ 
62.4 GeV. 
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TABLE IL Collision centrality selection in terms of percent- 
age of total Cu-Cu inelastic cross-section, number tracks in 
TPC, average number of participating nucleons and average 
number of binary nucleon-nucleon collisions for Cu-|-Cu at 
^/s^ = 200 GeV. 
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hadron multiplicity within the pseudo-rapidity range \r]\ 
< 0.5. In Table[T]we list the centrality bins for Au+Au at 
y^SNN — 62.4 GeV along with the multiplicity bin defini- 
tions, average number of participating nucleons and aver- 
age number of binary nucleon-nucleon collisions [ 29l . [30| . 
For the present analysis we chose six centrality bins cor- 
responding to 0-5%, 5-10%, 10-20%, 20-30%, 30-50%, 50- 
80% of the total inelastic nucleus-nucleus hadronic cross- 
section. A dataset of 2 million minimum-bias trigger 
events which passed the event cuts is used in the analysis. 
Tables HIl and IIIII list the six centrality bins for Cu-|-Cu 



TABLE III: Collision centrality selection in terms of percent- 
age of total Cu-Cu inelastic cross-section, number tracks in 
TPC, average number of participating nucleons and average 
number of binary nucleon-nucleon collisions for Cu+Cu at 
= 62.4 GeV. 
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at y^SNN = 200 GeV and 62.4 GeV corresponding to 0- 
10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60% of the 
total hadronic cross-section. The number of events used 
is 15 million and 24 niiUion for 62.4 GeV and 200 GeV 
Cu+Cu datasets, respectively, after the event cuts. 

C. Particle Selection 

We selected particle tracks in the rapidity region \y\ 
< 0.5. Particle identification was performed by corre- 
lating the specific ionization of particles in the TPC gas 
with their measured momenta. For this analysis pions 
are selected by requiring the specific ionization to be 
within 2 standard deviations from their theoretical Bich- 
sel value [Sll, [12] ■ In order to remove the kaons and pro- 
tons which could satisfy this condition, particles are also 
required to be more than 2 standard deviations from the 
Bethe-Bloch value for kaons and protons. Charged par- 
ticle tracks reconstructed and used for this analysis are 
accepted if they have space points on at least 15 pad rows 
in TPC. Tracks with fewer space points may be broken 
track fragments. These cuts are similar to those in our 
previous analysis of Au+Au collisions at ^snn — 200 
GeV [13] since the detector setup was identical. 

D. Pair Cuts 

Two types of particle track reconstruction errors di- 
rectly affect measured particle pair densities at the small 
relative momentum values studied here. Track splitting, 
in which one particle trajectory is reconstructed as two or 
more "particles," increases the apparent number of pairs 
at low relative q. To address this problem we developed 
a split track filter algorithm, described in our previous 
analysis of Au+Au collisions at ^snn — 200 GeV p7j . 



where values of the sphtting level parameter from —0.5 
to 0.6 ensured valid tracks. The inefficiencies arising 
due to track merging, in which two or more particle tra- 
jectories are reconstructed as one track, was completely 
eliminated by requiring that the fraction of merged hits 
(overlapping space-charge depositions in the TPC gas) be 
less than 10% for every track pair used in the correlation 
function. 

In the present analysis, we used the same cuts to re- 
move splitting and merging as were used for Au-f Au col- 
lisions at 



17[. The track pairs are 



= 200 GeV 

required to have an average transverse momentum (/ct 
= {\PiT + P2t|)/2) in one of 4 bins corresponding to 
[150,250] MeV/c, [250,350] MeV/c, [350,450] MeV/c and 
[450,600] MeV/c. The results are presented and discussed 
as a function of fcx and mx (= \/ fc^ 
those bins. 



mj) in each of 



as talai: 

C('?out, ^sidc, 9long) = (1^-^) + 
Ai^coul(*nv)(l+e-'?-*-«-* 



2 n2 2 r>2 

-^sidc-^sidc-Qlong-^lon 



0, (2) 



III. ANALYSIS METHOD 



where -ftTcoui is, to a good approximation, the squared 
nonsymmetrized Coulomb wave function integrated over 
a Gaussian source (corresponding to the LCMS Gaussian 
radii i?out, -Rsidc, ^long)- Assuming perfect experimen- 
tal particle identification and a purely chaotic (incoher- 
ent) source, lambda represents the fraction of correlated 
pairs f39| . 

We assumed a spherical Gaussian source of 5 fm for 
Au-|-Au collisions at ^snn = 62.4 GeV and a 3 fm source 
for Cu+Cu coUisions at ^/s^ = 62.4 and 200 GeV. The 
first term (1 - A) in Eq.(l2|) accounts for those pairs which 
do not interact or interfere and the second term repre- 
sents those pairs where both Bose-Einstein effects and 
Coulomb interactions are present [17] . 



Correlation function 



C. Systematic Uncertainties 



The numerator and denominator of the two particle 
correlation function in Eq.(IT]) are constructed by filling 
histograms corresponding to particle pairs from the same 
event and from mixed events, respectively. The back- 
ground pairs are constructed from mixed events where 
by pairing each particle in a given event is mixed with all 
particles from other events within a subset of ten simi- 
lar events. The events for mixing are selected within the 
given centrality bin such that their respective primary 
vertex z positions are all within 10 cm of one another. 



B. Bertsch-Pratt Parametrizations and Coulomb 
interactions 



We studied several sources of systematic errors simi- 
lar to a previously published STAR pion interferometry 
analysis for Au-|-Au collisions at ^snn = 200 GeV [17| . 
The following effects are considered: track merging, track 
splitting, source size assumed for the Coulomb correction, 
particle identification purity, and particle pair acceptance 
effects for unlike-sign charged pions. The estimated sys- 



tematic errors are less than 10% for i?niit, R' 



■side : 



R 



■long 1 



A 



in all centrality and /cx bins for the present datasets and 
are similar to those in p7| . This similarity is expected 
since the detector setup was identical and similar particle 
and pair selection cuts are used for Au-|-Au and Cu-fCu 
collisions. Results shown in the figures for the present 
datasets include statistical errors only. 



We decompose the relative momentum ~q accordin g to 
the Bertsch-Pratt (or "out-side-long" ) convention pHll^. 
[isl . [ssl . [3^ ] . The relative momentum ~q is decomposed 
into the variables (/long along the beam direction, gout 
parallel to the transverse momentum of the pair fcx = 
(pit + P2t)/2, and ^sido perpendicular to giong and gout- 

In addition to the correlation arising from the quantum 
statistics of two identical (boson) particles, correlations 
can also arise from two-particle final state interactions 
even for non-identical particles j55, 36, 37]. For identical 
pions the effects of strong interactions are negligible, but 
the long range Coulomb repulsion causes a suppression 
of the measured correlation function at small ~q . 

In this paper we follow the procedure used in our 
previous analysis of Au-|-Au collisions at ^/snn — 200 
GeV [13] . For an azimuthally integrated analysis at mid- 
rapidity in the longitudinal co-moving system (LCMS) 
the correlation function in Eq. ([1]) can be decomposed 



IV. RESULTS AND DISCUSSION 



A. Au+Au collisions at 



62.4 GeV 



The correlation function in Eq.© is fitted to the 
3D correlation data for Au-|-Au collisions at ^/skn = 
62 A GeV for each centrality and tot bins as defined 
above. The analysis is performed separately for 7r+7r+ 
and 7r~7r~ pairs. The final histograms for the like-sign 
pairs do not show appreciable differences and may there- 
fore be summed in order to increase statistics. Figure [T] 
presents the results for i?out, ^sidc, -Riong, A and the ratio, 
-Rout/^sido- The three femtoscopic radii increase with in- 
creasing centrality as expected, whereas the values of A 
and the i?out/^sido ratio exhibit no clear centrality de- 
pendences. 

We observe that for all centralities the three femto- 
scopic radii decrease with increasing tut whereas the A 
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FIG. 1: (Color Online) The femtoscopic parameters vs. mx 
for 6 different centralities for Au+Au collisions at ^snn = 
62.4 GeV. Only statistical errors are shown. The estimated 
systematic errors are less than 10% for -Rout, inside, i?iong, A 
in all centrality and bins. 
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FIG. 2: (Color Online) The comparison of femtoscopic mea- 
surements of Au+Au collisions at -^snn = 200 GeV and 62.4 
GeV for 0-5% most central events. Only statistical errors are 
shown for Au+Au collisions at y^SNN = 62.4 GeV. The es- 
timated systematic errors for Au+Au collisions at ^snn = 
62.4 GeV are less than 10% for -Rout, i?sidc, i^iong, A in 0-5% 
most central events and all fcx bins. The 200 GeV results are 
from [ll. 



FIG. 3: (Color Online) The energy dependence of femto- 
scopic parameters for AGS, SPS and RHIC from Refs. 

51, 52, si, m, H, In, 113. Energy dependences of 
pion femtoscopic parameters for central Au+Au, Pb+Pb and 
Pb+Au collisions are shown for mid-rapidity and (fcr) ~ 0.2- 
0.3 GeV/c. Error bars on NA44, NA49, CERES, PHOBOS 
and STAR results at ^snn = 130 and 200 GeV include sys- 
tematic uncertainties; error bars on other results are statisti- 
cal only. Only statistical errors are shown for Au+Au colli- 
sions at y^SNN = 62.4 GeV; the estimated systematic errors 
are less than 10% for Tiout, -Rside, -Riong. The PHOBOS results 
from [531 for ^snn = 62.4 and 200 GeV are slightly shifted 
horizontally for visual clarity. 



parameter increases with mx. Such behavior is consis- 
tent with our previous measurements at -^/snn = 200 
GeV • The increase of parameter A with ttit is due to 
the decreasing contribution of pions produced from long- 
lived resonance decays at higher transverse momenta. 
For comparison, in Fig. [2] we show the results for Au+Au 
collisions at -y/sNN = 62.4 GeV and 200 GeV for the most 
central collisions. We observe that the -Rout values are 
similar for both cases, but there are differences between 
the values of -Rsidc and -Riong- The Rout /Rside ratio de- 
creases with increasing mx , but the values are higher for 
= 62.4 GeV than for = 200 GeV. 

The observed dependences of the three femtoscopic 
radii are qualitatively consistent with models with col- 
lective flow [10, im |43|. Collective expansion results 
in position-momentum correlations in both transverse 
and longitudinal directions. In an expanding source 
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the correlation between the space-time points where 
the pions are emitted and their energy-momentum pro- 
duce a characteristic dependence of femtoscopic radii on 
TOT B m, [13, m, S, m, m, 113, Sg. The decrease in 
the "out" and "side" components can be described by 
models including transverse flow [13, 113, 113, SH], and 
the decrease in the "long" component by those with lon- 
gitudinal flow [13, S, SI, SI . 



B. Energy dependence of femtoscopic radii 

In Fig. [3] we present the energy dependences of the 
three femtoscopic radii and the ratio i?out/^sidc for the 
available data from AGS, SPS and RHIC. The results are 
compiled for Au-t-Au, Pb-|-Pb and Pb-|-Au collisions at 
mid-rapidity and for (fcx) ^ 0.2-0.3 GeV/c. The present 
measurements for Au-|-Au collisions at -^snn = 62.4 GeV 
are also included. The quality of the present STAR data 
with respect to statistical and systematic errors is signifi- 
cantly better than that reported by PHOBOS [lOl at the 
same energy. PHENIX results are not included because 
they were reported for broader centrality bins. WA97 
results are also omitted because they were measured at 
higher transverse momenta. 

Comparative studies are a necessary part of searches 
for nontrivial structures in the excitation function which 
might arise from a possible phase transition [3]. The ra- 
dius parameter i?sidc has the most direct correlation with 
the source geometry whereas i?out encodes both geometry 
and time scale information. Experimental results show 
that i?side decreases at AGS energies and then displays 
a modest rise with collision energy from SPS to RHIC. 
i?iong increases with collision energy after an initial de- 
crease at the lower AGS energies. For i?out the changes 
are very small. 

Hydrodynamic model calculations [3, [3| predict an 
enhancement in the ratio of -Rout/-Rsidc with increasing 
beam energy. The experimental results show no such be- 
havior. The measured ratios are better reproduced by 
the AMPT (A Multi-Phase Transport) model [H, how- 
ever the individual predicted radii have a steeper decrease 
compared to the experimental data [3]- An alternative 
model using a relativistic quantum mechanical treatment 
of opacity and the refractive index is capable of reproduc- 
ing the observed results |59|] , but strongly depends on the 
assumed initial conditions and neglects the time depen- 
dence of the corresponding optical potential. Hydrody- 
namic calculations [60| including viscosity offer another 
possible explanation for the above deviation between the 
data and model calculations as recently shown in [6l| . 
According to recent hydrodynamic calculations, the fem- 
toscopic radii can be described either by using the initial 
Gaussian density profile [131 or by including the combi- 
nation of several effects including: pre-thermal acceler- 
ation, a stiffer equation of state, and additional viscous 
corrections |63|. Other recent studies with a granular 
source model |64| also obtain a better description of the 
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FIG. 4: (Color Online) Femtoscopic parameters vs. ttit for six 
centralities for Cu-|-Cu collisions at ^snn ~ 62.4 GeV. Only 
statistical errors are shown. The estimated systematic errors 
are less than 10% for i?out, -Rsidc, ^long, ^ in all centrality and 
kT bins. 



experimental measurements of pion femtoscopic radii. 



C. Cu+Cu collisions at ^/s^ = 62.4 and 200 GeV 

The correlation functions are similarly constructed for 
Cu-hCu collisions at y/s^ = 62.4 GeV and 200 GeV. 
The extracted femtoscopic radii, -Rout, -Rside and i?iong, 
along with the A parameter and the ratio -Rout / -Rsido are 
presented in Figs. [H and E] for the 62.4 and 200 GeV 
data, respectively. The results are presented for six dif- 
ferent centralities and four tot bins. The highest bin 
[450 - 600] MeV/c of the most peripheral centrality (50 
- 60 %) in Cu-fCu collisions at ^snn = 62.4 GeV is 
omitted due to inadequate statistics for decomposition 
with the Bertsch-Pratt parametrization. For both col- 
lision energies the three femtoscopic radii increase with 
increasing centrality whereas the A parameter shows no 
centrality dependence. The tot dependences of the fem- 
toscopic radii are similar to that for Au-|-Au collisions. 
The -RoutZ-Rside ratios exhibit no clear centrality depen- 
dences for either energy. 



D. Comparison of femtoscopic radii for Cu+Cu 
and Au+Au collisions 

In Fig. [S] the femtoscopic source parameters A, -Rout, 
-Rsidc, -Riong and the ratio -Rout /-Rsido for central (0-5%) 
Au+Au collisions at -^snn — 200 GeV jT3| are compared 
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FIG. 5: (Color Online) Femtoscopic parameters vs. ttit for six 
centralities for Cu+Cu collisions at -^snn = 200 CeV. Only 
statistical errors are shown. The estimated systematic errors 
are less than 10% for i?out, ^?sidc, -Riong, ^ in all centrality and 
fcx bins. 



with central (0-10%) Cu+Cu collisions at same beam en- 
ergy. As expected, the femtoscopic radii for Cu+Cu col- 
lisions are smaller than for Au+Au collisions at the same 
beam energy. It is interesting that the values of the ratio 
^out/^sido for the two systems are similar. 

In Fig. [7] we extend the comparison of femtoscopic 
source parameters to include central (0-5%) Au+Au col- 
lisions at ysiiT = 62.4 GeV, central (0-10%) Cu+Cu col- 
lisions at y/suN — 62.4 and 200 GeV, and central (0-15%) 
TT+TT"'' and 7r~7r^ correlations from Au+Au collisions at 
62.4 GeV from the PHOBOS experiment [50]. The fem- 
toscopic radii for Cu+Cu coUisions at ^/suN — 62.4 GeV 
are smaller than those for Au+Au collisions at the same 
beam energy. The femtoscopic radii for Cu+Cu central 
collisions are similar for both energies. The variation of 
the i?out/-Rsidc ratio with tot is similar for the Au+Au 
and Cu+Cu collision data. 

In Fig. [5] we present the tot dependences of the ra- 
tios of femtoscopic radii for the most-central Au+Au and 
Cu+Cu collisions at ^/snn = 200 and 62.4 GeV. Ra- 
tios for the same colliding ion systems are close to unity 
whereas ratios of radii for Au+Au to Cu+Cu collisions 
are ~1.5. Although the individual radii decrease signifi- 
cantly with increasing mx the ratios in Fig. [5] show that 
the femtoscopic radii for Au+Au and Cu+Cu collisions 
at 62.4 and 200 GeV share a common tjit dependence. 
This result can be understood in terms of models [H, [6^ 
which use participant scaling to predict the femtoscopic 
radii in Cu+Cu collisions from the measured radii for 
Au+Au collisions at -^snn = 200 GeV, assuming the 



FIG. 6: (Color Online) The comparison of system size de- 
pendence in femtoscopic measurements of STAR Au+Au and 
Cu+Cu collisions at ^snn = 200 GeV. Only statistical errors 
are shown for Cu+Cu collisions at ^/swn = 200 GeV. The 
estimated systematic errors for Cu+Cu collisions at y^SNN = 
200 GeV are less than 10% for flout, Aside, Along, A in 0-10% 
most central events and fcx bins. The Au+Au results are 
from [l3l. 
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FIG. 7: (Color Online) The comparison of femtoscopic mea- 
surements of STAR Cu+Cu collisions at ^snn = 200 and 
62.4 GeV and Au+Au collisions at ^/s^ = 62.4 GeV. Only 
statistical errors are shown for STAR results. The estimated 
systematic errors for STAR results are less than 10% for -Rout , 
Aside, Along, A lu all ceutrality and fcx bins. The PHOBOS re- 
sults [sol for positive and negative pions in Au+Au collisions 
at y^SNN = 62.4 GeV are compared with STAR results. 
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FIG. 8: (Color Online) Ratios of femtoscopic radii at top 
centralities for Au+Au and Cu+Cu collisions at ^snn = 200 
and 62.4 GeV vs. niT- Only statistical errors are shown for 
Au+Au collisions at ^snn = 62.4 GeV and Cu+Cu collisions 
at -^snn = 62.4 and 200 GeV. The estimated systematic er- 
rors for Au+Au collisions at -^/snn = 62.4 GeV and Cu+Cu 
collisions at ^/swn ~ 62.4 and 200 GeV are less than 10% for 
RoMt, Rsidc, -Riong in all centrality and fcx bins. The 200 GeV 
results are from [17l |. 



radii are proportional to A-'^/^, where A is the atomic 
mass number of the colliding nuclei. 



E. Volume estimates and multiplicity scaling 

Estimates of the pion freeze-out volume Vf in terms 
of the femtoscopic radii are provided by the following 
expressions: 



Vf OC RontRsidcRl 



ong- 



(3a) 



(3b) 



However, the correlation lengths (femtoscopic radii) de- 
crease with increasing mx corresponding to an ttit de- 
pendent region of homogeneity which, in expanding 
source models, is smaller than the true collision volume 
at freeze-out. The volume estimates [Eqs. (|5a)) and (j3bp ] 
are obtained from the lowest mT bin, corresponding to 
the /ct region from 150 to 250 MeV/c as discussed in 
Sec. (IiTdI) . 

The Vf measurements using Eq. (|5a|) as a function of 
^snn are presented in Fig. [9] for Au+Au, Pb+Pb and 
Pb+Au collisions at mid-rapidity and for the lowest 
bin defined above. The results show two distinct do- 
mains: First, at the AGS where the volume measure de- 
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FIG. 9: (Color Online) The energy dependence of the pion 
freeze-out volume for heavy ion collision data from the AGS, 
SPS and RHIC estimated using Eq. pa|l . The references are 
given in the caption of Fig.|31 Only statistical errors are shown 
for Au+Au collisions at y/swN ~ 62.4 GeV. The estimated 
systematic errors for Au+Au collisions at ^snn ~ 62.4 GeV 
are less than 10% for i?out, -Rside, Rion^ in all centrality and 
kT bins. 
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FIG. 10: (Color Online) Pion freeze-out volume estimates 
as a function of number of participants and charged parti- 
cle multiplicity density for Au+Au at -^snn = 200 and 62.4 
GeV. Only statistical errors are shown for Au+Au collisions 
at y^SNN = 62.4 GeV. The estimated systematic errors for 
Au+Au collisions at ^snn = 62.4 GeV are less than 10% for 
-Rout, -Rsidc, i^iong iu all Centrality and At bins. The 200 GeV 
Au+Au collision results are from [ij . The lines in each panel 
represent linear fits to the data. 
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FIG. 11: (Color Online) Pion freeze-out volume estimates as 
a function of charged particle multiplicity density for Au+Au 
and Cu+Cu collisions. Only statistical errors are shown for 
Au+Au collisions at -^snn = 62.4 GeV and Cu+Cu collisions 
at y^SNN = 62.4 and 200 GeV. The estimated systematic er- 
rors for Au+Au collisions at ^snn = 62.4 GeV and Cu+Cu 
collisions at -^snn = 62.4 and 200 GeV are less than 10% for 
.Rout, -Rsidc, i^iong iu all ccutrality and A;t bins. The 200 GeV 
Au+Au collision results are from [13] . The lines in each panel 
represent linear fits to the data. 



creases, and second, in the SPS and RHIC energy regimes 
where a monotonic increase is observed. 

A detailed description of this non-trivial behavior was 
suggested in [6^ based on the hypothesis of constant 
mean free path length of pions at freeze-out. The expla- 
nation provided in [bG^ defines the pion mean free path 
length, A/, as: 



A, 



1 



(4) 



where p/ is the freeze-out density and a is the total 
cross-section for pions to interact with the surrounding 
medium. The freeze-out density can be expressed as the 
number of particles N in the estimated freeze-out vol- 
ume V/, divided by V/, resulting in the second expression 
in Eq. ([¥]). The denominator, TVct, can be expanded as 
the sum of the pion-pion and pion-nucleon contributions. 
At AGS energies the pion-nucleon term dominates since 
the pion-nucleon cross-section is larger than the pion- 
pion cross-section. Also, the number of nucleons at these 
lower energies at mid-rapidity exceeds the number of pi- 
ons. Hence, a decrease in the number of mid-rapidity 
nucleons leads to a decrease in the observed freeze-out 
volume iyj) as a function of ^snn- At SPS and RHIC 
energies the pion-pion term dominates the denominator 
in Eq. ([4]) due to copious pion production leading to an 
increase in the observed Vf. 

Based on this interpretation we expect the volume es- 
timates in the pion dominated RHIC regime to show a 
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FIG. 12: (Color Online) The pion source radii dependences on 
charged particle multiplicity density for Au+Au and Cu+Cu 
collisions. Only statistical errors are shown for Au+Au col- 
lisions at y^SNN = 62.4 GeV and Cu+Cu collisions at -^snn 
= 62.4 and 200 GeV. The estimated systematic errors for 
Au+Au collisions at y^SNN = 62.4 GeV and Cu+Cu collisions 
at -ysNN = 62.4 and 200 GeV are less than 10% for -Rout, RsiAa, 
-Riong in all centrality and fcT bins. The 200 GeV Au+Au col- 
lision results are from [13]- The lines represent linear fits to 
the data. 



linear dependence on charged particle multiplicity. In 
Fig. [TO] freeze-out volume estimates (using Eqs. ([5a|) and 
(|3bp ) are shown as a function of the number of par- 
ticipants (left panels) and charged particle multiplicity 
(right panels) for Au+Au collisions at ^snn — 62.4 and 
200 GeV. The predicted linear increase with charged par- 
ticle multiplicity is observed. Estimated freeze-out vol- 
umes for Au+Au collisions at the same centralities in- 
crease with collision energy indicating that A^part is not a 
suitable scaling variable in this case. On the other hand, 
charged particle multiplicity provides better scaling prop- 
erties. 

Additional estimates of freeze-out volume dependences 
on charged particle multiplicity are presented in Fig. [11] 
for both the Au+Au and Cu+Cu results at -^snn = 62.4 
and 200 GeV. Both freeze-out volume estimates for the 
four collision systems show an approximate, common lin- 
ear dependence on charged particle multiplicity. The lin- 
ear dependences of femtoscopic radii on (dTVch/rfr;)^/'^ for 
Au+Au and Cu+Cu collisions at ^Snn = 62.4 and 200 
GeV are shown in Fig. 1121 The above common, linear 
dependences Q are consistent with the assumption of a 
universal pion mean-free-path length at freeze-out [66j . 
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V. SUMMARY AND CONCLUSIONS 

We have presented systematic measurements of pion 
femtoscopy for Au+Au collisions at ^snn = 62.4 GeV 
and Cu+Cu collisions at ^snn = 62.4 and 200 GeV, 
and compared these new results with our previous analy- 
sis of Au+Au coUisions at ^snn = 200 GeV [13 • For 
all the systems considered the three femtoscopic radii 
(i?out, -Rsidc and i?iong) increase with centrality, whereas 
the values of the A parameter and ratio i?out/^sidc are 
approximately constant with centrality. The three fem- 
toscopic radii decrease with increasing mx, whereas the 
A parameter increases with mx. The increase of A with 
Tot is attributed to decreasing contamination from pi- 
ous produced from long-lived resonance decays at higher 
transverse momentum. 

The decrease of femtoscopic radii with increasing tot 
can be described by models with collective, transverse 
and longitudinal expansion or flow. The ratios of femto- 
scopic radii at top centralities for different colliding sys- 
tems (Au-I-Au and Cu-t-Cu) at ^snn = 62.4 and 200 
GeV show that the corresponding radii vary similarly 
with Tot. 

The predicted rise of the ratio i?out/-Rsidc with col- 
lision energy due to a possible phase transition 3] is 
not observed for Au-|-Au and Cu-|-Cu collisions. The 
compilation of freeze-out volume estimates Vf as a func- 
tion of collision energy ^/snn (using Eq. (|5a)) along with 
the datasets presented in Fig. [3]) shows two distinct do- 
mains: with increasing ^snn, Vf decreases at the AGS, 
but steadily increases throughout the SPS and RHIC en- 
ergy regimes. At AGS energies the decreasing number 
of baryons at mid-rapidity leads to a decrease in the ob- 
served freeze-out volume {Vf) as a function of y/suN- At 
higher beam energies from SPS to RHIC copious and in- 
creasing pion production causes the freeze-out volume to 



rise. 

The dependences of the freeze-out volume estimate on 
number of participants and charged particle multiplicity 
are compared. Measurements for Au-|-Au collisions at 
the same centralities, but different energies yield differ- 
ent freeze-out volumes demonstrating that A^part is not 
a suitable scaling variable. The freeze-out volume esti- 
mates for all four collision systems presented here show 
a linear dependence on final charged particle multiplic- 
ity which is consistent with the hypothesis of a universal 
mean- free-path length at freeze-out. 

For the systems studied here the multiplicity and kj^ 
dependences of the femtoscopic radii are consistent with 
previously established trends at RHIC and at lower ener- 
gies. The radii scale with the final state collision multi- 
plicity which, in a static model, is consistent with an hy- 
pothesized universal mean- free-path length at freeze-out. 
This and similar studies establish the baseline systemat- 
ics against which to compare future femtoscopic studies 
at the LHG ;67]. 
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